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ABSTRACT 

We present a deep, wide-field near-infrared (NIR) survey over five widely separated fields at high 
Galactic latitude covering a total of ~ 3 deg 2 in J, H, and K s . The deepest areas of the data (~ 
0.25 deg 2 ) extend to a 5 a limiting magnitude of JHK S > 24 in the AB magnitude system. Although 
depth and area vary from field to field, the overall depth and large area of this dataset make it one 
of the deepest wide-field NIR imaging surveys to date. This paper discusses the observations, data 
reduction, and bright galaxy counts in these fields. We compare the slope of the bright galaxy counts 
with the Two Micron All Sky Survey (2MASS) and other counts from the literature and explore the 
relationship between slope and supergalactic latitude. The slope near the supergalactic equator is sub- 
Euclidean on average pointing to the possibility of a decreasing average space density of galaxies by 
~ 10 — 15% over scales of ~ 250 — 350 Mpc. On the contrary, the slope at high supergalactic latitudes 
is strongly super-Euclidean on average suggesting an increase in the space density of galaxies as one 
moves from the voids just above and below the supergalactic plane out to distances of ~ 250 — 350 
Mpc. These results suggest that local large scale structure could be responsible for large discrepancies 
in the measured slope between different studies in the past. In addition, the local universe away from 
the supergalactic plane appears to be underdense by ~ 25 — 100% relative to the space densities of a 
few hundred megaparsecs distant. 

Subject headings: cosmology: observations and large scale structure of universe — galaxies: funda- 
mental parameters (counts) — infrared: galaxies 



1. INTRODUCTION 

Historically, deep and wide-area surveys in the near- 
infrared (NIR) have been mutually exclusive. Deep sur- 
veys covered only small solid angles while wide surveys 
covered large solid angles, to only bright limiting magni- 
tudes. Both types of surveys suffer from a limited range 
of apparent magnitudes over which good statistics could 
be obtained. The boundaries between the two have come 
closer together in recent years with the development of 
large format NIR imagers capable of covering large areas 
to significant depths. Bridging the gap between "deep" 
and "wide" NIR photometry is now becoming possible, 
allowing a coherent observational picture of the history 
of the universe and large scale structure to be developed. 

In this paper we present a dataset composed of five 
large and widely separated fields at high galactic lati- 
tudes, including one centered on an Abell cluster (Abell 
370). Here we present the observations, data reduction, 
and bright galaxy counts in these fields. In total, our 
survey consists of ~ 2.75 deg 2 to a 5 a limit in a 3" aper- 
ture of JHK S ^ 22 — 23 and an additional ~ 0.25 deg 2 
to JHK S ~ 24, making it one of the deepest wide-field 
surveys to date in the NIR. 

The NIR samples a region of the spectral energy dis- 
tribution (SED) of galaxies that is insensitive to galaxy 
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type over a large redshift range (|Cowie et~aT][l99l . As 
such, galaxy counts in the NIR can be used as a tracer 
of large scale structure without having to make assump- 
tions about the distribution of galaxy types in the counts. 
Bright galaxy counts are expected to be dominated by 
M < M* galaxies at low redshift (z < 0.2). The logarith- 
mic derivative ( "slope" hereafter) of the bright counts is 
expected to follow the Euclidean prediction (a = 0.6) 
until the knee of the luminosity function induces a tran- 
sition to a lowe r slope around an apparent magnitude of 
JHK S ~ 17.5 (jBarro et alJl2009f ). Thus, if one assumes 
no recent evolution in the NIR galaxy luminosity func- 
tion, then the slope of the counts curve is a measure of 
homogeneity in the local universe on scales of hundreds 
of megaparsecs (M* ~ 16 at 500 Mpc). 

Several authors have found evidence for a large 
(hundreds of megaparsecs in radiu s) local underden- 
sity of galaxies using galaxy counts ( Huang et al.l 119971 : 



sity o l galaxies using galax y counts (Huang e t al.lll»»rt 
iBusswell et ail [2001: iFrith etaTI |200j,l2QQ5|). An under- 
density on these scales could imply that local measure- 
ments of the Hubble consta nt are high by some tens of 
percent dHuang et al.lll997h. The Tw o Micron All Sky 
Survey (2MASS. iSkrutskie et al1l2006h probes to depths 
of JHK S ~ 15 providing a measure of very nearby large 
scale structure. A challenge to deeper NIR surveys has 
been to cover a large enough area of sky to achieve good 
counting statistics at magnitudes just beyond the reach 
of 2MASS. The wide-field aspect of our data allow us to 
investigate the question of a possible local void by linking 
our galaxy counts to those of 2MASS and calculating the 
slope of the galaxy counts curve in the local Universe. 

Catalogs for this work will be presented in an upcoming 
paper (R. Keenan et al. 2010, in preparation), where we 
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select one of the largest samples to date of Distant Re d 
Galaxies (DRGs, J - K Wcga . < 2.3, iFranx et al.ll200l . 
We will also use this NIR survey in conjunction with 
deeper K a — band data to constrain the NIR cosmic in- 
frared background. In addition, we are engaged in a 
campaign of spectroscopic followup of H— band selected 
galaxies from this survey to constrain the local galaxy 
mass function and merger rates to z ~ 1. Unless other- 
wise noted, all magnitudes given in this paper are in the 
AB magnitude system (t/iab = — 2.5 log./), — 48.60 with 
/„ in units of ergs cm~ 2 s _1 Hz -1 . 

The structure of this paper is as follows: In Section [2] 
and Section [3] we describe the fields observed, the expo- 
sure times, the limiting magnitudes, and the data reduc- 
tion methods. In Section [4] we describe the photometry 
methods and the completeness of the data. We present 
star-galaxy separation methods and galaxy counts in Sec- 
tion [5] In Section [6] we analyze our counts with those 
from 2MASS and other published studies in a discussion 
of galaxy counts and large scale structure. In Section [7] 
we summarize our results. 

2. OBSERVATIONS 

Our five NIR fields (CLASXS, CLANS, CDF-N, A370, 
and SSA13) overlap with existing data at other wave- 
lengths and cover a large and diverse area of sky, which 
minimizes the effects of cosmic variance and large scale 
structure. We give a summary of the exposure times, 
field areas, and limiting magnitudes in Table [T] Below 
are descriptions of our five fields and existing ancillary 
data for each. 

Our first two fields are centered on t he Chandra Large 
Area Synoptic X-ra y Survey (CLASXS; lYang et al.|[2004 : 
iSteffen et al.l [2004) an d the Chandra Lockman Area 
North Survey (CLANS: ITrouille et al.l[27M[2009h . Each 
of these fields cover ~ 1 dcg 2 in JHK S . These fields 
are located in the Lockman Hol e region of extremel y 
low Galactic HI column density ((Lockman ct al. 1986). 
CLASXS consists of nine overlapped ~ 40 ks Chandra 
exposures (with the central field observed for 73 ks) com- 
bined to create an ~ 0.4 deg 2 image, and CLANS con- 
sists of nine separate ~ 70 ks Chandra e xposures com- 
bined to cr eate an ~ 0.6 deg 2 image. iTrouille et al.l 
(|2008l 120091 ) present spectroscopic redshifts for over half 
of the > 1000 X-ray sources in these fields, and they 
present photometric redshifts for almost all of the re- 
maining sources. They also present optical data on 
these fields from CFHT MegaCam (u,g',i' for CLASXS 
and g',r',i',z' for CLANS) and Subaru SuprimeCam 
(B,V,R,I,z' for CLASXS). These fields are covered in 
the Spitz er Wide- Area Infrar ed Extragalactic Survey 
(SWIRE. [Lonsdale et al"1l2003f ) Legacy Science Program 
at 3.6 fim. and 24 /*m to limiting fluxes of 5 /uJy and 
230 /^Jy, respectively. A large portion of the CLANS field 
was also recently covered in a deep 20 cm Very Large Ar- 
ray (VLA) survey containing ~ 2000 sourc es to a limiting 
flux of ~ 20 /xJy jO wen fe Morrison|[2008f) . 

Our third field covers a 0.25 deg 2 area cen- 
tered on the Chandra Deep Field North (CDF-N, a 
0.12 deg 2 2 Ms Chandra observation; iBrandt et al.l 
[200lt lAlexander et~al] [20031 ). The entire NIR field 
has also been observed in the optical at B, V, R, I, 
and z' with Subaru SuprimeCam and in the [/—band 



with the KPNO 4 m MOSAI C ins trument (iCapak et al l 
20041) . iBarger et all (|2003f ) and ITrouille et all (|2008l . 



20091) present a highly complete spectroscopic redshift 



catalog for X-ray sources in the CDF-N along with corre 
sponding optical, NIR, and Spitzer Legacy Science Pro- 
gram 3.6 and 24 /im photometry (M. Dickinson et al. 
2009, in preparation). The CDF-N contains the Great 
Observatories Origins Deep Survey North (GOODS-N; 
145 arcmi n 2 HST Advance Ca mera for Surveys ob- 
servation, iGiavalisco et all 12004) with data at F345W, 
F606W, F775W, and F850LP. Barger et al. (2008, and 
references therein) present a highly complete spectro- 
scopic catalog of sources in the GOODS-N and deep 
K s — band photometry from CFHT WIRCam. In addi- 
tion, we have ultradeep K s — band image of the GOODS- 
N taken with the MOIRCS instrument on the Subaru 
Telescope. The GOO DS-N has also b een imaged at 
20 cm with the VLA bv lRichardsl (|2000D . iBiggs fc Ivisonl 
(2006), and G. Morrison et al. (2009, in preparation). 

Our fourth field is the Abell 370 (A370) cluster and 
surrounding area (~ 0.5 deg 2 ). A3 70 is a cluster of 
richness at a redshift of z = 0.37. We have unpub- 
lished optical data o n this field fr o m Su baru Suprime- 
Cam (B,V,R,I,z'). IBarger et all (|2001D have studied 
this field previously in optical, NIR, and radio. We have 
followed up spectroscopically thousands of sources in and 
around the A370 cluster, some of which have been pub 



lished for high-redshift Lyman alpha emitters (Hu et al 



2002) and for low-metallicity galaxies (jKakazu et al 



I2007D . The deep 20 cm VLA data and more of the spec- 
troscopic follow-up on this field will be presented in I. 
Wold et al. (2010, in preparation). 

Our fifth field is a ~ 0.2 deg 2 area centered on the 
"Small-Survey- Are a 13"(SSA13) from the Hawaii Deep 
Fields desc ribed inlLillv et all (119911) . iMushotzkv et al.l 
( 2000) and IBarger et al.l "( 20011 ) have observed Chandra 
X-ray sources in this field in optical, NIR, and submil- 
limeter wavelength s. The entire NIR field has also been 
imaged at 20 cm bv lFomalont et akl (|2006[). for which op- 
tical s pectr oscopy have been pr esented in lChapman et al.l 
(120031) and lCowie et all ((200l . 

2.1. ULBCam Observations 

We made all the H and J— band observations (except 
the CDF-N J— band) with the Ultra Low Background 
Camera (ULBCam) on the UH 2.2 m telescope. The 
ULBCam NIR imager was the first wide-field infrared 
camera. It combined four Hawaii2-RG HgCdTe detectors 
to cover a 17' x 17' field of view (FOV) with a plate scale 
of 0'.'25 per pixel. ULBCam was developed as a prototype 
for the NIR imager on th e James Webb Space Telescope 
(Jl^gr. lHaTl et al.ll200l . 

Our ULBCam observations began in September 2003. 
Observations have continued since then totalling over 300 
hours of observing through February 2008. All the data 
to date have been reduced and are mosaicked into the 
final images. Average seeing in the final mosaics is ~ 1". 

We used a 13-step dither sequence consisting of 45" — 
90" shifts between images having exposure times of 
30 — 120 s depending on the background levels. The 
maximum exposure time of 120 s kept the entire dither 
sequence less than 30 min such that the coadded images 
from one sequence would represent a relatively constant 
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TABLE 1 

Coordinates, exposure times, areas, and 5 a limiting magnitudes for each field 



Field CDF-N A370 CLANS CLASXS SSA13 

R.A.(hh:mm:ss) 12:36:55 02:39:53 10:46:54 10:34:58 13:12:16 

Dec (dd:mm:ss) 62:14:19 -01:34:58 59:08:26 57:52:22 42:41:24 

Galactic 1 (deg) 125.9 173.0 148.2 151.5 109.1 

Galactic b (deg) 54.8 -53.5 51.4 51.0 73.8 

Supergal. 1 (deg) 54.7 302.3 52.2 52.3 75.3 

Supergal. b (deg) 11.7 -25.7 -1.8 -3.9 13.6 



J— band observations 



Instrument CFHT WIRCAM UH 2.2m ULBCAM UH 2.2m ULBCAM UH 2.2m ULBCAM UH 2.2m ULBCAM 
Max. exposure time (hr) 7.7 11.2 4.8 4.9 5.7 
Area in square deg. 0.24 0.55 1.05 0.71 0.22 
3" aper. 5ct limit 24.4 23.5 23.1 23.0 23.2 



//—band observations 



Instrument UH 2.2m ULBCAM UH 2.2m ULBCAM UH 2.2m ULBCAM UH 2.2m ULBCAM UH 2.2m ULBCAM 
Max. exposure time (hr) 24.5 5.9 3.9 2.4 5.6 
Area in square deg. 0.4 0.59 1.24 1.1 0.27 
3" aper. 5ct limit 24.0 23.0 22.5 22.5 22.9 



K s — band observations 



Instrument CFHT WIRCAM UKIRT WFCAM UKIRT WFCAM UKIRT WFCAM UKIRT WFCAM 
Max. exposure time (hr) 36.8 3.1 15.1 5.9 4.0 
Area in square deg. 0.29 0.87 0.91 1.0 0.87 
3" aper. 5a limit 24.8 23.0 23.4 23.5 23.1 



sky background. 

2.2. WIRCam Observations 

We observed the CDF-N in the K s — band with the 
Widefield Infrared Camera (WIRCam) at the Canada 
France Hawaii Telescope 3.6 m (CFHT). We observed a 
large area (~ 0.25 deg 2 ) centered on the field. WIRCam 
combines four HAWAII2-RG 2kx2k detectors that cover 
a 20' x 20' FOV with a pixel scale of 0'.'3. We dithered 
images to cover the detector gaps and to provide a uni- 
form sensitivity distribution. We performed most of the 
observations under photometric conditions with seeing 
from 0'.'6 — 1". The final mosaic seeing is ~ 0'.'9. 

Our group made the K s — band observations in 
semesters 2006A and 2007A, and we combined these 
data with public observations of the GOODS-N made 
with the same instrument by a Canadian group led by 
Luc Simard in 2006A. Altogether the final K s image in- 
cludes ~ 40 hours of integration time, making this one 
of the deepest wide-field K s -band images ever (5 a limit 
K s ~ 24.8). The jfg-b and catalog for this field is given 
in iBarger et all (|200l) . 

J— band observations with WIRCam were obtained by 
a group led by Lihwai Lin in 2006A. Our group reduced 
these public data in 2008. The final mosaic includes ~ 9 
hours of integration time with mosaic seeing ~ 0'.'8 and 
a 5 a limiting magnitude of J ~ 24 . 

2.3. WFCam Observations 

We observed four of our five fields in K— band with 
WFCam on UKIRT (the CDF-N is inaccessible to 



UKIRT because of declination constraints) . WFCam em- 
ploys four HAWAII2-RG 2kx2k detectors spaced by 94%. 
This design, combined with a pixel scale of 0'.'4, allows 
for coverage of ~ 0.77 deg 2 in a mosaic of four pointings. 
Average seeing in the final WFCam mosaics is ~ 1 — 1'.'2. 

3. DATA REDUCTION 

Although the details of our data reduction varied de- 
pending on the instrument used, the overall process was 
essentially the same. High sky backgrounds with strong 
variability in the NIR require that many relatively short 
exposures be coadded to create deep images. Dithered 
pointings on high Galactic latitude fields ensure that the 
sky background can be effectively modelled and sub- 
tracted from mosaicked frames. Many short exposures 
also allow for the removal of cosmic rays by determining 
where transient signals appear from frame to frame. 

We reduced all the ULBCam and WIRCam data using 
the Interactive Data Language (IDL) SIMPLE Imaging 
and Mosaicking Pipeline (SIMPLE; W.-H. Wang 2009, in 
preparation). This pipeline was designed by W.-H. Wang 
specifically for use in reducing dithered wide-field near- 
infrared images from ground-based mosaic cameras. In 
its current form, the pipeline is optimized for blank-field 
extragalactic surveys where there are no large extended 
objects. Below we provide a synopsis of the pipeline pro- 
cedures, but for extensive documentation and the code 
itself we refer the reader to W.-H. Wang's website 5 . 

In the case of WFCam, we retrieved partially reduced, 

5 http:/ /www.asiaa.sinica.edu.tw/~whwang/idl/SIMPLE 
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stacked images from the WFCam Science Archive (WSA) 
and applied the final steps of the SIMPLE pipeline to do 
additional background subtraction, as well to determine 
the absolute astrometry and photometry for each stacked 
frame before combining them into a final mosaic. 

For calibration we adopted the 2MASS zero magnitude 
fluxes, which are 1594 Jy, 1024 Jy and 666.7 Jy for J, H, 
and K s , respectively. The conversion from Vega to AB 
for 2MASS magnitudes are Jab = Jvcga + 0.894, Hab = 
Hy esa + 1.374, and if Sj AB = if s ,vcga+ 1-84. We analyzed 
all of our final images for flatness by looking at a flux 
ratio between our data and 2MASS for point sources as 
a function of position on each image. In this analysis we 
found that our final mosaics are quite flat and uniform 
with no discernable systematic differences in the flux ra- 
tio as a function of position of the image and with typical 
random errors of < 0.05 magnitudes. 

In this paper we use the catalogs f rom the United 
King dom Infrared Sky Survey (TJKIDSS. lLawrence et al.l 
12007 ) for calibration checks and galaxy counts compari- 
son. UKIDSS uses the UKIRT Wide Field Camera (WF- 
CAM, |C^al^t^Lll2003L^he photometric system is de- 
scribed in Hcw ett et al.1 (120061). and the calibration is de- 
scribed in Hodgkin et al.l (j2009f l. The pipeline processing 
and science archive are described in M. Irw in et al. (2009, 
in preparation) and lHamblv et a l. ( 2008). We have used 
data from the 3rd release, which will be described in de- 
tail in S. Warren et al. (2009, in preparation). 

3.1. ULBCam Reduction Details 

We reduced images within each dither set (typically 15- 
30 minutes in length) using the SIMPLE pipeline. SIM- 
PLE features a robust method for flat fielding in which 
a sky flat is iteratively derived from dithered night sky 
images. Application of this method leads to extremely 
flat images after background subtraction. SIMPLE also 
corrects for image distortion in a set of dithered images 
without any prior knowledge about the optics and with- 
out the use of an astrometric catalog, which allows for 
accurate registration of wide-field images. 

We applied the following procedures within the SIM- 
PLE pipeline to generate our reduced mosaics: First, 
we derived a median sky flat and appli ed it to all frames. 
We th en used the SExtractor package (|Bertin fc Arnoutsl 
1996) to detect objects in each flattened frame. We 
masked detected objects to create a new median sky flat 
and then redid the flatfielding. We subtracted the sky 
background by fitting a smooth polynomial surface to 
each masked, flattened frame. We removed the brightest 
cosmic ray hits with a 5x5 pixel spatial sigma filter. 

Next, we again used the SExtractor package to mea- 
sure object positions and fluxes in each flattened, sky- 
subtracted image. We calculated the first-order deriva- 
tive of the optical distortion function by measuring the 
offsets of each object in the dither sequence as a function 
of location in the image. We obtained absolute astrome- 
try by matching detected objects to a reference catalog. 
For this purpose we used t he catalogs from t he Sloan Dig- 
ital Sky Survey (SDSS) (jYork et al.l l2000h for all fields 
but Abell 370, where USNO B1.0 was used due to the 
lack of SDSS coverage. We then warped the reduced ex- 
posures to a common tangential sky plane with subpixel 
accuracy. This projection corrects for both optical dis- 
tortion and absolute astrometry. Through this process 



we achieved rms astrometry error of ~ O'.'l — Of.' 2, except 
in the Abell 370 field where with the USNO B1.0 catalog 
we achieve an rms error of ~ 0'.'5. 

We then weighted all projected images by their relative 
SExtractor object fluxes to correct for variable extinction 
and combined them to form a dither mosaic. In this 
combination we applied a sigma filter to pixels that have 
the same sky position to further remove faint cosmic rays 
and artifacts such as reflections inside the optics. 

We applied SExtractor a third time to each dither mo- 
saic and scaled the object fluxes into /xJys (resulting in an 
AB zcropoint of 23.9 fo r each image) th rough a compari- 
son with 2MASS (|Skrutskie et al.ll2006f) for point source 
photometry in each field in the range 14 < JH < 16. 
This is the range over which 2MASS reports better than 
S/N — 10 a for point sources and our data are well be- 
low saturation levels. In the final step we coadded all 
reduced dither mosaics into a large final mosaic. 

3.2. WIRCam Reduction Details 

The reduction procedure for the WIRCam data (which 
is only the CDF-N field) was the same as for the ULB- 
Cam data, except for the following details: WIRCam 
shows significant crosstalk in raw frames between the 32 
readout channels. This was removed by subtracting the 
median of the 32 64 x 2k channels in the object-masked 
image. This process removed the majority of crosstalk, 
though weak features persist around the few brightest 
stars in the field. 

We obtained absolute astrometry for the WIRCam 
frames by matching detected objects to a reference cat- 
alog constructed with relat ively bright and compa ct ob- 
jects in the ACS catalog (Giavalisco et al. 2004; after 
correcting for the C/4 offset between AC S and radio 
frame s) and the SuprimeCam catalog of Capa k et alJ 
(2004). The detection limits of 2MASS in the if s -band 
only marginally overlap with the linear detection regime 
of WIRCam, and so rather than using 2MASS, we cali- 
brated WIRCam K s — band fluxes in each individual re- 
duced fram e using our deeper S ubaru MOIRCS image as 
a reference (jBarger et al.ll2008l ). In the J-band 2MASS 
is deeper and we were able to achieve a good calibration 
for WIRCam frames using 2MASS point sources in the 
range 14 < J < 16. 

Although the WIRCam and ULBCam reduction tech- 
niques were similar, the data come from different cameras 
on different telescopes. As a check for consistency, we 
compared the final reduced CDF-N J— band WIRCam 
image with the CDF-N J— band mosaic we generated 
from the ULBCam data. The result was that for 62 rel- 
atively bright stars (14 < J < 16) we found a magnitude 
offset of Jwir — Jvlb — 0.003 ±0.05, and for 691 galaxies 
at J < 20 we found J wm - J ULB = 0.05 ± 0.05. These 
small systematics are of the same order of the random 
errors in the data, and so we consider reduced data from 
these two cameras to be consistent with one another. 

3.3. WFCam Data Reduction 

We elected to use the stacked WFCam frames from the 
Cambridge Astronomical Survey Unit (CASU) pipeline, 
rather than the micro-dithered raw frames, to take ad- 
vantage of the existing WFCam reduction pipeline and 
the "dribbling" algorithm employed by CASU to bring 
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the relatively undersampled (0'.'4 per pixel) raw frames 
to higher effective resolution (C'2 per pixel). The fol- 
lowing is a brief overview of the WFCam data reduction 
pipeline. Extensive documentation of the pipeline can be 
found at the CASU- WFCam data processing website 6 . 

WFCam data are preprocessed in the UKIRT summit 
pipeline, in which raw frames are corrected for linearity 
and instrument signature, then shipped to CASU for fur- 
ther processing. The CASU pipeline includes astromet- 
ric alignment, sky background subtraction, and stacking 
of micro-stepped images. First, dark frames are stacked 
to generate a master dark which is subtracted from all 
frames. Next, a flatfielding is performed using a stacked 
master twilight flat. A sky frame is then generated us- 
ing a median combination with sigma clipping of a se- 
ries of dark sky science frames. This is an iterative pro- 
cess including object detection and masking to form the 
most robust sky background estimation. The resulting 
sky background is then subtracted from science frames. 
Individual exposures are both dithered (tens of arcsec- 
onds offsets) and micro-stepped (subpixel offsets). The 
micro-stepping allows for the final step in the pipeline of 
"dribbling" (similar to drizzling) to stack microstepped 
raw frames to a finer grid. Processed frames are archived 
at the WSA. 

We retrieved the stacked frames from the WSA and 
applied the post-stacking steps of the SIMPLE pipeline 
described above to perform additional background sub- 
traction and spatial sigma filtering to get rid of any resid- 
uals from the initial reduction. We obtained absolute as- 
trometry by matching detected objects in each stacked 
frame to a reference catalog. We then warped the re- 
duced stacks to a common tangential sky plane with sub- 
pixel accuracy to correct for both optical distortion and 
absolute astrometry. Through this process, we achieved 
rms astrometry error of ~ O'.'l — 0'.'2, except in the Abell 
370 field where with the USNO B1.0 catalog we achieve 
an rms error of ~ 0'.'5. 

Our CLANS field overlaps with the UKIDSS Deep Ex- 
tragalactic Survey (DXS), and as such we retrieved the 
public DR3 data from the WSA and added it to our 
final CLANS if-band mosaic. The UKIDSS DXS ob- 
servations employed a very similar approach to our ob- 
servations and could be combined directly with our data 
using the same reduction techniques. 

Finally, we applied SExtractor to each dither mosaic 
and scaled the object fluxes into /LtJys (resulting in an 
AB zeropoint of 2 3.9 for each image) th rough a compari- 
son with 2MASS (|Skrutskie et al.ll2006T) for point source 
photometry in each field in the range 14 < K s < 16. In 
the final step we coadded all reduced dither mosaics into 
a large final mosaic. 

Regarding the calibration of our WFCam data using 
2MASS, the WFCam K filter is a classical K filter with 
a central wavelength of ~ 2.2 /im, while the 2MASS K s 
filter has a central wavelength of ~ 2.16 fim. The dif- 
ference between photometry for stars in these two fil- 
ters can be empirically quantified, as is done in the gen- 
eration of UKIDSS catalogs, where 2MASS K s is used 
to calibrate WFCam K data via the following relation: 

if WFCam = ^s,2MASS + 0.01 ( J2MASS — ^s,2MASs)- 

We use the above UKIDSS empirical relation to cali- 

6 http:/ /casu. ast.cam.ac.uk/surveys-projects/wfcam/technical/ 



brate our WFCam data using 2MASS. However, for the 
analysis in this paper we wish to compare our WFCam 
K— band galaxy counts with 2MASS K s — band counts. 
The K s — K magnitude difference for galaxies is differ- 
ent than that of stars and depends both on the shape 
of the galaxy's SED and its redshift. To explore the 
possible systematic and random errors associated with 
generating K s — band galaxy counts from K— band data, 
we employed two m ethods. 

First, we used the lBruzual & CharlotJ (|2003f ) template 
SEDs for a wide range of stellar population ages. We 
find that out to a redshift of ~ 2.5, the magnitude differ- 
ence K s — K has a median value of ~ —0.1 with rms 
variability of ±0.05. The difference is slightly larger 
(~ — 0.15±0.05) at the lowest redshifts. Second, we com- 
pared magnitudes from our own photometry of bright 
galaxies in our four WFCam fields (described in Sec- 
tion |4|) to 2MASS magnitudes for the same objects. We 

find a K s -K difference of 0.15±0.1 between 2MASS 

and our magnitudes, in goo d agreement with the low - 
rcdshift estimation using the iBruzual fc Chariot! ([2003ft 
templates described above. Thus, in our galaxy counts 
and analysis, we adopt a zeropoint offset to our if— band 
data of —0.15 magnitudes to bring the counts in line with 
the 2MASS if s -band. 

In the case of our WFCam if —band fields, we did not 
have photometry from ULBCam or WIRCam to compare 
with directly as a check on the different reduction meth- 
ods. We did, however, check our photometry against the 
UKIDSS catalogs in the CLANS field where we overlap 
and found a magnitude offset for 564 bright (15 < K < 
18) stars and galaxies of -Kukidss — K = 0.03 ± 0.05 
which, again, is of the same order as the random errors 
in the data. In addition, the fact that our if —band pho- 
tometry when compared with 2MASS K s is in agree- 
ment with the prediction s from simulations using the 
IBruzual &; Chariot! (|2003[ ) templates suggests the reduc- 
tions have performed as expected. 

4. PHOTOMETRY 

We generated catalogs for all fields in our three 
NIR bands using SExt ractor software version 2.4.4 
(jBertin fc Arnouta 119960 for source identification and 
photometry. Parameters used in the SExtrac- 
tor configuration file are given in Table [2] Pa- 
rameters listed as "variable" in the "setting" col- 
umn (DETECT.MINAREA and SEEING.FWHM) were 
calculated individually for each field. The DE- 
TECT.MINAREA parameter sets the minimum num- 
ber of contiguous pixels required to be considered a 
valid source. We reduced DETECTJVUNAREA until 
a spike in the faint magnitude galaxy counts was ob- 
served (due to spurious detections) and set the final value 
just above this limit. The SEEING_FWHM parameter 
is used by SExtractor to morphologically classify sources 
as pointlike or extended. We calculated our input SEE- 
ING.FWHM values as described in Section O DE- 
TECT.MINAREA and SEEING_FWHM values for each 
field are listed in Table El MAG.ZEROPOINT values 
were set to 23.9 in all cases corresponding to units of 
/zJy in the images. 

We masked areas near bright objects to exclude areas 
of possible saturation or diffraction spikes. We define 
bright objects to have a USNO-B1.0 R— band magni- 
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TABLE 2 
S Extractor Parameters 



Parameter 


Setting a 


Comment 13 


TiT7 ir ri?r~ ir r tvdp 
Dcj 1 HiKj J Li rUt 


L A_. \ ) 


Detector type 


DH; 1 rL/O 1 _IV111N ArtrAjA 


Variable 


Mininmum number of pixels above threshold 


TATTTTrr^T TUD T7CTU 


1 


Detection threshold (<r) 


A T\T A T VQTQ TTTRPQU 
A IN i\Li i 1 JrLri jzjoxI 


i 

i 


Limit for isophotal analysis (cr) 


T7TT TT7D 


v 


Use filtering (Y or N) 


TTTT T'TTR M A A/TT7 
r lijl JzjJ\._1n AIVijij 


gauss_2.0_5x5.conv 


Filter for detection 


■p\T7"D T TPT\T"Pl TvT r PTU"D T?QU 

JJriir>L;r J lNiJ_iN Ixlxti^brl 


32 


Number of deblending sub-thresholds 


Dri/BLrLlNU JVUJNCUJN 1 


0.001 


Minimum contrast parameter for deblending 




Y 


Clean spurious detections (Y or N) 


CLrL AIM _rAHAM 


1.0 


Cleaning efficiency 


AT A QTV TVDT? 




Correct flux for blended objects 


DTJ/^T A T)T7 , TZ> r PT TTD T7C 

.FHO 1 _AFrLHl UKrLo 


3, 6 


MAG_APER aperture diameters in arcseconds 


r nU 1 _AU 1 UirVVxlAJVlO 


z.o, o.o 


MAG_AUTO parameters: Kron_factor,min_radius 


C A TT T"D T T?"\7"I7T 


OUUUU.U 


Saturation level 


HAP 7T7T5/~\Di^MM r P 

MALj_z>ii/KtJFUllN 1 


23.9 


Magnitude zeropoint 


(~* A TM 


1.0 


Gain is 1 for absolute rms map 


PIXELJ3CALE 





size of pixels (0=use FITS WCS info) 


SEEING_FWHM 


Variable 


Stellar FWHM in arcseconds 


STARNNW_NAME 


default, nnw 


Neural- Network weight table filename 


BACK.SIZE 


48 


Background mesh in pixels 


BACK_FILTERSIZE 


8 


Background filter 


BACKPHOTO.TYPE 


LOCAL 


Photometry background subtraction type 


WEIGHT.GAIN 


N 


Gain does not vary with changes in rms noise 


WEIGHT.TYPE 


MAP.WEIGHT 


Weighting with exposure map 


MEMORY.OBJSTACK 


2000 


number of objects in stack 


MEMORY_PIXSTACK 


200000 


number of pixels in stack 


MEMORY_BUFSIZE 


1024 


number of lines in buffer 



a Settings listed as "variable" were calculated individually for each field in each band. Values for the 
DETECT.MINARE and SEEING.FWHM are given for each field in all bands in Tabic [3] 
k The SExtractor manual containing details of the SExtractor parameter function is available through 
TERAPIX at http://tcrapix.iap.fr 



TABLE 3 

Input SEEING_FWHM and DETECT.MINAREA Parameters to SExtractor for Each Field in JHK s 



Field J H K 



SEEING_FWHM a DETECT_MINAREA a SEEING.FWHM DETECT.MINAREA SEEING.FWHM DETECT.MINAREA 
CDF-N 0.88 9 1.20 9 0.92 4 
A370 1.20 12 1.28 9 1.08 12 
CLANS 1.20 12 1.16 16 1.12 9 
CLASXS 1.12 9 1.24 16 1.20 9 
SSA13 1.16 16 1.20 12 1.20 16 



a Calculated as described in Section^ 

5. GALAXY COUNTS 

Our survey covers ~ 3 deg 2 in J and H and nearly 
4 deg 2 in K s . However, only 1.8 deg 2 is covered simul- 
taneously in all 3 bands. As shown in Section 15.11 we 
reliably separate stars from galaxies at JHK S < 17.5 by 
morphology alone. As such, for JHK S < 17.5 we count 
galaxies over the entire extent of our survey to achieve 
the best counting statistics possible. For JHK S > 17.5 
we switch to counting galaxies over only the 1.8 deg 2 of 
overlap in all three bands in order to also employ color- 
selection techniques in separating stars from galaxies. 

We constructed the galaxy counts by first making sev- 
eral cuts to the output catalog of SExtractor. We ex- 
cluded objects in masked areas as described in Section|4j 
We identified stars by the methods described in Sec- 
tion 15.11 The raw counts for JHK S < 24 are shown 
in Figure [I] and the raw counts for JHK S < 22, where 
all fields are highly complete, are given in tabular form 
in Tables [ffl 



tude < 14.5. Following the method used in lCapak et al.l 
pool , we mask an area around such objects within a 
radius defined by (where R stands for R— band magni- 
tude): 

Radius = 291.2 - 34.55 x R + 1.06 x R 2 (1) 

We also masked regions of low exposure time in the mo- 
saic. We defined low exposure time to be any area having 
less than 50% of the median exposure time in the mosaic. 

As a measure of the depth of our fields, we calculated 
5 a limiting magnitudes (listed in Table [l} by measuring 
the rms noise in randomly scattered 3" apertures at blank 
locations in the image. We did this by scattering 10 4 
apertures randomly and retrieving rms noise values for 
those that fell at least 3" away from detected objects. We 
retrieved the median rms value in all scattered apertures 
and recorded this as the noise background of the image 
(i.e., 1 cr). 
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Fig. 1. — (a-c) Raw galaxy counts in each of our five fields. J— band counts are shown in (a), H— band in (b) and K s — band in (c). Error 
bars represent 1 <r Poisson fluctuations, (d-f) The same data as in (a-c) divided through by an arbitrary normalized Euclidean model to 
expand the ordinate. 



We binned identified galaxies in variable sized bins us- 
ing large (1 magnitude) bins at bright magnitudes and 
decreasing to small (0.25 magnitude) bins at faint mag- 
nitudes. This scheme of decreasing the bin size toward 
fainter magnitudes allowed us to provide better statis- 
tics on the bright end of the galaxy counts curve and 
better resolution on the faint end. We divided the bin to- 
tals by the unmasked area of the image. We determined 
Poissonian galaxy c ount errors (1 a) using the counting 
confidence limits of Gehrelsj (|1986f ). 

We then performed a weighted average of galaxy counts 
over all fields. We created the weighted average by tak- 
ing into account Poisson errors (weight = l/o- 2 ). This 
ensures that results from the widest fields in our survey 
carry more weight on the bright end. 

5.1. Star-Galaxy Separation 

We used a combination of morphology and color to 
separate stars from galaxies for our galaxy counts. Our 
methods include use of the SExtractor "CLASS_STAR" 
parameter, the ratio of semimajor to semiminor axes of 
objects, and J — K s color to separate stars from galaxies. 
We honed our methods on a group of spectroscopically 
identified stars and galaxies in the GOODS-N (a subfield 
of the CDF-N). This involved the use of J and K s — band 
catalogs fr om this study and t he GOODS-N redshift cat- 
alog from Ba rger et alj (|2008D . Finally, as an indepen- 
dent test of our methods, in the four of our five fields (all 



but A370) where we overlap with the SDSS, we compared 
our classified stars having JHK S < 19 with the morpho- 
logical classifications for these objects in the SDSS cata- 
logs. The results of these methods are discussed below. 
Our final star counts generated using these methods are 
shown in Figured 

The CLASS_STAR parameter is an output of SExtrac- 
tor to describe whether a given object looks like a point 
source or an extended source. The CLASS_STAR out- 
put contains a value between (extended object) and 1 
(point source) and is assigned by SExtractor to each ob- 
ject detected in the image. CLASS_STAR can be used 
to separate stars from galaxies at bright magnitudes, but 
the accuracy of this parameter breaks down at fainter 
magnitudes. Where that breakdown occurs depends on 
the depth of the field and the SEEING_FWHM input 
parameter (see Table [5]). 

The SEEING_FWHM input parameter is provided 
to SExtractor by the user to describe the average 
point source Full- Width-Half-Maximum (FWHM) in 
the image. This parameter is critical for SExtractor 
to accurately identify stars with its output parameter 
CLASS_STAR. 

We determined the optimal value for the SEE- 
ING_FWHM input parameter by first running SExtrac- 
tor on an image to determine the FWHM of point sources 
in the image (using an initial guess at the image seeing 
for the input SEEING.FWHM). We then took the me- 



Keenan et al. 



TABLE 4 
J — band Raw Galaxy Counts 



Mag(AB) log(JV) mag~ 1 deg- 2 ± lo- a 



CDF-N CLANS CLASXS SSA13 A370 



15. 


.000 


1 


,134± 


0. 


361 





,105± 


0. 


.518 


0. 


236± 


0. 


518 




0± 0. 


265 


1 


133± 


0. 


225 


16 


.000 


1. 


310± 


0. 


293 


1 


.717± 





.063 


1. 


651± 


0. 


077 


1 


,751± 





.163 


1 


,756± 


0. 


085 


17. 


.000 


2 


,230± 


0. 


079 


2 


303± 





.033 


2 


,241± 


0. 


041 


2 


260± 


0. 


.073 


2 


,483± 


0. 


039 


17. 


.750 


2 


,564± 


0. 


067 


2 


,744± 





.035 


2. 


,716± 


0. 


035 


2 


698± 


0. 


.060 


2 


,954± 


0. 


031 


18. 


.250 


2. 


,907± 


0. 


046 


2 


939± 


0. 


.028 


2. 


862± 


0. 


030 


2 


,974± 





.045 


3 


,203± 


0. 


023 


18. 


.750 


3 


,054± 


0. 


039 


3 


,178± 





.022 


3 


,092± 


0. 


023 


3 


,119± 


0. 


.038 


3 


,371± 


0. 


019 


19. 


.250 


3. 


321± 


0. 


029 


3 


391± 


0. 


.017 


3. 


339± 


0. 


.017 


3 


,270± 


0. 


.032 


3. 


,506± 


0. 


016 


19. 


.750 


3 


593± 


0. 


021 


3 


,568± 


0. 


.014 


3 


,550± 


0. 


014 


3 


,522± 


0. 


.024 


3 


.701± 


0. 


013 


20 


.250 


3. 


,763± 


0. 


018 


3 


,800± 


0. 


.010 


3. 


.750± 


0. 


Oil 


3 


,708± 





.020 


3. 


,851± 


0. 


Oil 


20 


.750 


3. 


,962± 


0. 


014 


3 


,962± 


0. 


.009 


3. 


971± 


0. 


008 


3 


,925± 


0. 


.015 


4 


033± 


0. 


009 


21 


.125 


4. 


060± 


0. 


018 


4 


109± 


0. 


.010 


4. 


115± 


0. 


010 


4 


,077± 


0. 


.018 


4. 


,183± 


0. 


010 


21 


.375 


1 


159± 


0. 


016 


4 


,200± 


0. 


.009 


1 


218± 


0. 


009 


4 


,208± 





.015 


1 


,273± 


0. 


009 


21 


.625 


1 


,228± 


0. 


015 


4 


,277± 


0. 


.008 


4. 


.344± 


0. 


008 


4 


,374± 


0. 


.013 


1 


,363± 


0. 


008 


21 


.875 


4. 


311± 


0. 


013 


4 


379± 


0. 


.007 


4 


,448± 


0. 


007 


1 


436± 





.012 


4 


,444± 


0. 


008 



a Errors represent Poissonian confidence limits of lCekrell ll^bl) 



TABLE 5 

H — band Raw Galaxy Counts 



Mag(AB) log(A0 mag- 1 deg~ 2 ± lo- a 



CDF-N CLANS CLASXS SSA13 A370 



15. 


.000 


0. 


,609± 


0. 


518 





,983± 


0. 


.163 


1. 


.071± 


0. 


163 




0± 0. 


265 


1 


,038± 


0. 


255 


16. 


.000 


1 


,563± 


0. 


163 


1 


,792± 





.053 


1. 


,719± 


0. 


063 


1 


,948± 





.124 


1 


834± 


0. 


079 


17. 


.000 


2. 


,281± 


0. 


059 


2 


399± 


0. 


.027 


2. 


.290± 


0. 


034 


2 


395± 





.062 


2. 


,566± 


0. 


035 


17. 


.750 


2. 


,813± 


0. 


051 


2 


,818± 


0. 


.033 


2. 


.725± 


0. 


035 


2 


,868± 


0. 


.050 


3. 


031± 


0. 


028 


18. 


.250 


3 


,034± 


0. 


040 


3 


065± 





.025 


2. 


,925± 


0. 


028 


3 


,042± 


0. 


.041 


3 


,261± 


0. 


022 


18. 


.750 


3. 


,278± 


0. 


031 


3 


,262± 


0. 


.020 


3. 


133± 


0. 


022 


3 


,169± 





.036 


3 


,380± 


0. 


019 


19. 


.250 


3. 


,520± 


0. 


023 


3 


,458± 


0. 


.016 


3. 


375± 


0. 


017 


3 


419± 


0. 


.027 


3 


,555± 


0. 


015 


19. 


.750 


3 


696± 


0. 


019 


3 


665± 





.012 


3 


,576± 


0. 


013 


3 


613± 


0. 


.022 


3 


,750± 


0. 


012 


20 


.250 


3. 


,917± 


0. 


015 


3 


,882± 


0. 


.009 


3. 


.797± 


0. 


010 


3 


,835± 





.017 


3. 


,907± 


0. 


010 


20 


.750 


1 


,080± 


0. 


012 


4 


,051± 


0. 


.008 


3. 


.986± 


0. 


008 


i 


,011± 


0. 


.014 


4 


,138± 


0. 


008 


21 


.125 


4. 


193± 


0. 


015 


4 


,187± 





.009 


4. 


131± 


0. 


010 


4 


,191± 


0. 


.016 


4. 


300± 


0. 


009 


21 


.375 


1 


,289± 


0. 


014 


4 


,281± 


0. 


.008 


4. 


.224± 


0. 


009 


i 


,298± 





.014 


4 


,397± 


0. 


008 


21 


.625 


4. 


375± 


0. 


012 


4 


366± 





.008 


4. 


,312± 


0. 


008 


4 


365± 





.013 


1 


,492± 


0. 


007 


21 


.875 


4. 


,444± 


0. 


011 


4 


,443± 





.007 


4. 


,394± 


0. 


007 


4 


419± 


0. 


.012 


4 


,585± 


0. 


006 



a Errors represent Poissonian confidence limits of lGehrefsl 1198611 



TABLE 6 

K — band Raw Galaxy Counts 



Mag(AB) log(A0 mag~ 1 deg~' 2 ± lo- a 



CDF-N CLANS CLASXS SSA13 A370 



15. 


.000 


1. 


,047± 


0. 


361 


1 


,028± 


0. 


.176 


0. 


.976± 


0. 


176 





900± 





.204 


1 


139± 


0. 


155 


16 


.000 


1 


649± 


0. 


176 


1 


,951± 





.050 


1. 


,763± 


0. 


057 


1 


,794± 


0. 


.059 


2. 


003± 


0. 


048 


17. 


.000 


2 


,444± 


0. 


057 


2 


,634± 





.023 


2. 


,585± 


0. 


023 


2 


,575± 


0. 


.025 


2 


752± 


0. 


020 


17. 


.750 


2. 


,934± 


0. 


045 


2 


,914± 


0. 


.029 


2. 


818± 


0. 


032 


2 


,955± 


0. 


.046 


3. 


,182± 


0. 


024 


18. 


.250 


3 


,188± 


0. 


034 


3 


,152± 





.022 


3 


,038± 


0. 


025 


3 


,107± 


0. 


.039 


3 


375± 


0. 


019 


18. 


.750 


3 


,425± 


0. 


026 


3 


,418± 





.016 


3 


,317± 


0. 


018 


3 


396± 


0. 


.028 


3. 


.502± 


0. 


016 


19. 


.250 


3. 


,683± 


0. 


019 


3 


609± 


0. 


.013 


3. 


.507± 


0. 


014 


3 


602± 


0. 


.022 


3. 


,698± 


0. 


013 


19. 


.750 


3 


,864± 


0. 


016 


3 


,821± 


0. 


.010 


3 


,722± 


0. 


011 


3 


736± 


0. 


.019 


3 


,846± 


0. 


011 


20 


.250 


1 


,018± 


0. 


013 


3 


,978± 


0. 


.008 


3. 


890± 


0. 


009 


3 


932± 





.015 


3. 


,995± 


0. 


009 


20 


.750 


4. 


110± 


0. 


012 


4 


108± 





.007 


4. 


,021± 


0. 


008 


4 


,087± 


0. 


.013 


4. 


,147± 


0. 


008 


21 


.125 


4. 


252± 


0. 


014 


1 


219± 





.009 


4. 


130± 


0. 


010 


4 


239± 


0. 


.015 


4. 


,243± 


0. 


010 


21 


.375 


4 


,271± 


0. 


014 


4 


,288± 


0. 


.008 


4. 


222± 


0. 


009 


4 


,276± 


0. 


.014 


1 


,329± 


0. 


009 


21 


.625 


4 


,370± 


0. 


012 


4 


,344± 


0. 


.008 


4. 


.298± 


0. 


008 


4 


,385± 





.013 


1 


,372± 


0. 


008 


21 


.875 


4. 


423± 


0. 


012 


4 


,412± 





.007 


4. 


,378± 


0. 


007 


4 


,420± 


0. 


.012 


4. 


,449± 


0. 


008 



a Errors represent Poissonian confidence limits of |Oehrell |1§53) 
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Fig. 2. — Star counts for each of our fields using the com bination 
of morphology and color criteria described in Section 15.11 




16 18 20 22 24 
K s (AB) 

Fig. 3. — J — K s color-magnitude plot for spectroscopically iden- 
tified stars and galaxies in t he GOODS-N region (a subfield of the 
CDF-N) from the catalog of Bargcr ct al. (2008). Accurately iden- 
tified stars are shown in blue. Stars we fail to identify are shown 
in purple. Accurately identified galaxies are shown in red circles 
and galaxies misidentified as stars are in green circles. The vast 
majority of stars are found below the line J — K a = 0.1. 




-1.0 -0.5 0.0 0.5 1.0 1.5 
J-K s 

Fig. 4. — A color-color plot demonstrating the utility of a combi- 
nation of the SDSS g' filter with the J and K s filters in the NIR for 
star-galaxy separation. Stars are shown as bl ue stars and galaxie s 
as red circles from the spectroscopic catalog of Bargcr ct al. (2008). 
We use the g' JK S color-color data from our fields that overlap with 
the SDSS to investigate discrepancies between our sta r-gal axy sep- 
aration and that of the SDSS, as described in Section |5. II 



dian FWHM for selected point sources (CLASS.STAR 

> 0.95) as the final input value for the SEEING.FWHM 
parameter. Values calculated in this way for the SEE- 
ING-FWHM parameter in each field are listed in Table 

m 

We found that on occasion for bright galaxies the 
CLASS_STAR parameter misidentified galaxies as stars 
that were clearly galaxies by eye (2 out of 125 galaxies at 
K s < 19 in the GOODS-N catalog). We resolved this is- 
sue by constraining the ratio of scmimajor to scmiminor 
axis for stars (a/b < 1.3). 

Using the star identification criteria of CLASS-STAR 

> 0.7 for JK S < 17.5 and CLASS-STAR > 0.9 and a/b < 
1.3 for 17.5 < JK S < 19, we correctly identified 100% of 
stars for JK H < 1 9 from the spectroscopic sample of 
IBarger et all (f2008h . 

To extend our star-galaxy separation to fainter mag- 
nitudes, we used a combination of the SExtractor 
CLASS-STAR parameter and the J — K s color of ob- 
jects. We determined that a color cut of J — K s < 0.1 
and CLASS-STAR > 0.95 for 19 < JK S < 24 accu- 
rately identified > 90% of stars while only misclassifying 
< 1% of galaxies as stars. Because galaxies dominate 
the counts at these faint magnitudes, the combination of 
these effects introduces less than 1% error in the counts 
in the range 19 < JK S < 24. These results are dis- 
played in Figure [3] where K s — band magn i tude i s plotted 
against J — K s color for the Bargcr ct al. (2008) sample. 
Blue stars show correctly identified stars, red circles show 
correctly identified galaxies, and purple stars and green 
circles show misclassified stars and galaxies, respectively. 

In our fields that overlap with SDSS (all but A370), we 
were able to perform a check on our bright (JHK S < 19) 
star selection by cross-correlating with objects classified 
morphologically as point sources in SDSS. We looked at 
two things in this analysis. First, we searched for ob- 
jects that we identifed as stars which were classified as 
extended objects in SDSS. Second, we searched for ob- 
jects that we found to be galaxies that were classified as 
point sources in SDSS. 

For JHK S < 17.5, we found that none of our galaxies 
were listed as point sources in SDSS. However, a handful 
(5-15) of our stars in each field were listed as extended 
sources in SDSS. We first investigated this discrepancy 
by eye and found that many of these objects were clearly 
point sources in our images identifiable by their diffrac- 
tion spikes. This is perhaps due to the fact that our 
images are much deeper than SDSS, or these are late- 
type stars that are brighter in the NIR than in the op- 
tical. Other sources appeared to be close binary point 
sources in our images, perhaps unresolved by the SDSS. 
Given that the majority of discrepant point sources ap- 
pear to be misclassifications in the SDSS, we conclude 
that we have done a robust job of star-galaxy separation 
for JHK S < 17.5. 

In the range 17.5 < JHK S < 19, we found > 90% 
agreement between objects we classfied as stars and point 
sources identified by SDSS, as well as for objects we clas- 
sified as galaxies and their SDSS counterparts. To in- 
vestigate discrepant sources, we employed g'(sdss) —J 
versus J — K s color-color diagrams. The utility of such 
diagrams for sepa rating stars and ga laxies is shown in 
Figure [4] using the IBarger et all (|2008fl sample as an ex- 
ample, where SDSS counterparts are found for 131 stars 



10 



Keenan et al. 



and 647 galaxies. Using this sample, we determined that 
a color cut of g' — J > 6 . 5 x ( J — K s ) + 1 . 8 was appropriate 
for selecting the vast majority of stars. 

We plotted objects we identified as stars, but whose 
SDSS counterparts were classified as extended on the 
g'JK s plane and found that the majority matched the 
color criterion for stars defined above. Those that did 
not meet the criterion typically had very high SExtractor 
CLASS_STAR parameters in our catalogs and, as such, 
are likely candidates to be late- type M and L dwarf stars. 

Finally, we plotted objects that we identified as galax- 
ies but whose SDSS counterparts were classified as point 
sources on the g'JK s plane. We found a minority of 
these objects did, in fact, meet the g'JK s color criterion 
for stars. In this case, we flagged these objects as stars in 
our catalogs. Stars identified in this way through com- 
parison with the SDSS were very few in number and only 
affect the galaxy counts in the 17.5 < JHK S < 19 range 
on the ~ — 3% level. 

The SDSS does not cover our A370 field, so we do not 
have the opportunity to compare our star-galaxy sepa- 
ration with an independent optical morphology classifi- 
cation like we have done for the other fields. Based on 
our comparison with the SDSS in other fields, we expect 
that stars we may have missed in the A370 field could 
contribute to an excess on the few percent level to counts 
in the 17.5 < JHK S < 19. range. 

As seen in Figure [3j there are stars that will be missed 
by our color cut in J — K s . These will tend to be M and 
L dwarf stars because they are redder than the rest of 
the main sequence population (T dwarfs, in fact, have 
J — K s < like typical main sequence stars). However, 
we would expect to identify the majority of M and L 
dwarfs as stars based on morphology alone for K s < 19 
(as demonstrated above), corresponding to distances of 
out to ~ 100 parsecs. Beyond this we would miss these 
faint dwarf stars in our star removal. 

The space density of M and L dwarfs is poorly 
constrained a t larg e he liocentric dis t ances , but 
iStanwav etHI (|2008t ) and iPirzkal et all (j2009f ) have 
shown that M and L dwarfs may exist in large numbers 
out to several kpc into the Galactic halo. The star 
counts for M dwarfs from these two groups roughly 
agree and appear flat out to very faint magnitudes at 
levels of - 100 - 200 mag" 1 deg 2 . 

In our galaxy counts, assuming the worst case scenario, 
this could represent a 15% systematic error in the bright- 
est bin (K s = 19 — 20), for which we require J—K s < 0.1. 
Assuming flat star counts for these late-type dwarf stars, 
this error would drop to < 10% at K s — 20.5 and < 5% 
at K s — 21.5. The focus of this paper is the bright 
galaxy counts for which we believe our star-galaxy sepa- 
ration is robust, but we alert the reader that the counts 
remain uncorrected for late-type dwarf stars fainter than 
JHK S = 19. 

6. BRIGHT GALAXY COUNTS AND LOCAL LARGE SCALE 
STRUCTURE 

In the past two decades, numerous NIR galaxy counts 
studies have been published. Typically new data are 
overplotted with a selection of previous works to demon- 
strate how the various studies agree or disagree in a given 
bandpass. At bright magnitudes, factors of two or more 
discrepancy between any two studies are not uncommon 
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Fig. 5. — NIR K— corrections are shown as a function of red- 
shift i n J (black), H (red), and K s (green), from Mannucci et al. 
120011 '). Five galaxy types (E, SO, Sa, Sb, Sc) are shown (though not 
identified individually) to demonstrate the similarity in each NIR 
bandpass of the K— correction for different galaxy types. Because 
of this similarity the uncertainty in type mixture in our galaxy 
counts only leads to uncertainties of < 0.005 in the determination 
of the expected Euclidean slope. 
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Fig. 6. — The distance of galaxies at a given K 3 — band mag- 
nitude (the distributions are very similar in H and J— band) as- 
suming homogeneity and a non-evolving NIR luminosity function. 
This distance dist ribution was calcu lated using the NIR lumi- 
nosity functions of I Jones et al.l (2006) and the K— corrections of 
IMannucci et al] 120011) . The plotted points show the peak of the 
distance distribution (M*) in a given magnitude bin, and the er- 
ror bars show one standard deviation above and below the peak in 
each bin. 



in these comparisons. A lack of knowledge about the 
proper normalization of the galaxy counts curve makes it 
difficult to determine when such differences are actually 
due to variations in large scale structure and not other 
effects, such as differences or inconsistencies in photom- 
etry. 

The slope of the bright galaxy counts curve, on the 
other hand, is a measure of local large scale structure 
that is independent of the normalization. In a static, 
homogeneous, Euclidean universe, galaxy counts would 
follow the curve J\f = 4x 10 Qm , where A is a normal- 
ization constant, a = 0.6, and m is apparent magnitude. 
A slope of a = 0.6 is also expected for the bright end of 
the galaxy counts curve in an expanding, homogeneous 
universe where K— corrections are small and nearly in- 
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dependent of galaxy type and where cosmological effects 
o n the luminosity di stance are small. 

iBarro et al.l (120091) s how that for galaxy counts derived 
from iSchechterl (|1976[ ) luminosity functions, the bright 
counts slope is expected to follow the Euclidean predic- 
tion and be dominated by M < M* galaxies until a 
transition to lower slope is induced by the knee in the 
luminosity function around JHK S ~ 17.5. 

The bright counts slope is then modified by the 
if— corrections in a given bandpass, which act effectively 
as a changing M*. Optical if —corrections are large 
and galaxy type dependent, and thus considerations of 
large scale structure based on the slope of optical galaxy 
counts are plagued with uncertainties about the exact 
mix of galaxy types in the counts. NIR if— corrections 
are smaller than in the optical and nearly type inde- 
pend ent. Figure [5] shows the low-redshift if— corrections 
from rMannucc i et al.l ((2001) (J in black, H in red, and 
K s in green). These are shown for five galaxy types (E, 
SO, Sa, Sb, Sc), demonstrating the near type indepen- 
dence of the NIR if —corrections. 

It is straightfo rward to determine t he effect of the 
if —corrections of iMannucci et aD (|2001h on the theoret- 
ical slope of the galaxy counts curve in a homogeneous 
distribution of galaxies in the nearby expanding universe 
(z < 0.1, cz = H Q D). Assuming H = 70 km s _1 Mpc -1 
and an equal mix of all five galaxy types shown in Fig- 
ure [SJ the slope for the J— band should be a ~ 0.57, for 
the H— band a ~ 0.58, and for the K s — band, a ~ 0.61. 
Because of the similarity of the if —corrections across 
galaxy types the uncertainty introduced by the lack of 
knowledge of the exact type mixture in the calculation 
of the expected slope is less than ±0.005. This result is 
also independent of the shape of the galaxy luminosity 
function if one assumes no galaxy evolution at z < 0.1. In 
Figure [6] we show the distribution in distance of galaxies 
for a range of if s — band magnitudes. 

In what follows, we consider departures of ~ 0.01 — 0.1 
from the Euclidean prediction . If we ignore evolution 
and treat the slope as a direct measure of local large scale 
structure, these departures can be thought of as changes 
in the space density of galaxies as one moves further out 
into the local universe. For example, a difference of 0.01 
from the Euclidean predication measured over an inter- 
val of three magnitudes would roughly imply a galaxy 
space density change of ~ 7% over the range of distances 
sampled at those magnitudes. As shown in Figure [6l the 
range of distances sampled in any given magnitude bin 
increases toward fainter magnitude, adding uncertainty 
to the physical significance of a change in slope, but if 
the luminosity function is truly a constant over the range 
measured, then the slope is a measure of the space den- 
si ty of galaxies oyer th at range. 

iMaddox et al.1 (|1990l ) found a steeper than Euclidean 
slope in B— band counts over the 4300 deg 2 of the Auto- 
mated Plate Measuring (APM) galaxy survey and sug- 
gested rapi d galaxy evolution at low redshifts as an ex- 
planation. IKoo fc Kronl (|1992l ) point out that dramatic 
collective evolution near the present epoch would be un- 
usual, and they instead attribute the steep slope of the 
APM survey counts to some co mbination of m odel uncer- 
tainties and normal evolution. iShanksl (|1990l ) point out 
that the steep galaxy counts slope could be due to clus- 



tering if we are located inside a large (~ 150/i _1 Mpc) 
region of ~ 50% underdensity. However, as mentioned 
above, optical if — corretions are highly type dependent, 
so a measurement of the slope in the NIR is required for 
a more reliable test o f local large scale structure. 

iHuang et al.l (|1997f ) found a steeper than Euclidean 
slope (a — 0.69) in if s — band counts over a collection 
of fields totalling ~ 10 deg 2 , but they point out that 
the amount of evolution needed to account for the ef- 
fect, particularly in the NIR at low redshifts, seems 
unlikely. The alte rnative to rapid evolution cited by 
IHuang et al.1 (|l997p, is a ~ 50% local underdensity of 
galaxies stretching 300ft. -1 Mpc in extent. They point 
out that this would mean local measurements of the cos- 
mological mass density, Oo, would be low by the same 
factor and that local measurements of Hq could be high 
by as much as 33%. Evidence for such a "local void" us- 
ing 2MASS, the Two-degree Field Galaxy Redshift Sur- 
vey (2dFGRS), and other dat a, along with ga l axy co unts 
models, has been presen ted in iBusswell et al.l (|2004l ) and 
iFrith et al.l (|2003l . 120051 ). In what follows, we look at the 
bright galaxy counts slope of our data in combination 
with 2MASS, the variability in the slope as a function of 
position on the sky, and possible implications for local 
large scale structure. 

6.1. Galaxy Counts in the 2MASS XSC 

We downloaded the entire 2MASS extended source cat- 
alog (XSC) and 2MASS-6x XSC from the NASA/IPAC 
Infrared Science Archive 7 . We constructed 2MASS 
galaxy counts from the catalog over an area of ~ 20, 000 
deg 2 representing the entire 2MASS XSC for Galac- 
tic latitude \b\ > 30 (2ir Steradians). We limited this 
2MASS sample to JHK S < 15.5 to stay within the all 
sky catalog's completeness limits. We excluded known 
2MASS Galactic extended sources from the counts us- 
ing the Galactic extended source catalog available at 
NASA/IPAC website. As shown in Figure [6j galaxy 
counts from 2MASS from 11-15*' 1 magnitude are sam- 
pling M* galaxies in the local universe at distances of 
~ 50 - 350 Mpc. 

We also constructed ga laxy counts from the 2MASS- 
6x (|Beichman et al.ll2003l ) catalog of extended sources in 
the ~ 25 deg 2 centered on t he Lockman Hole are a of low 
galactic HI column density (|Lockman et al.lfl986T ) ■ These 
data represent 25 deg 2 where 2MASS observed six times 
as long as the all sky catalog exposure times, thus mak- 
ing the 2MASS-6x catalog roughly one magnitude deeper 
in J,H and if s . We were able to further "deepen" the 
2MASS-6x Lockman Hole catalog because our deep data 
in the CLANS and CLASXS fields overlaps 100% with 
the 2MASS-6x Lockman Hole. As such, we were able 
to measure completeness in the 2MASS-6x catalog, and 
correct the counts to JHK S = 16.75. Thus, we use the 
2MASS and 2MASS-6x catalogs to firmly establish aver- 
age galaxy counts on the bright end. The average 2MASS 
galaxy counts for |6| > 30 deg are shown in Figure [71 

The solid lines in FigureUJshow an error-weighted least 
squares fit to the data. In the left panels the counts are 
shown in the standard N vs. magnitude form, and in the 
right panels the counts are shown divided through by a 

7 http://irsa.ipac.caltech.edu 
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Fig. 7. — (a-c) Average counts for 2MASS |b| > 30 (black diamonds) and for 2MASS-6x Lockman Hole (blue triangles). Error bars 
show 1 a Poisson fluctuations (smaller than plot symbols in a-c). (d-f) are the same data as (a-c) after having been divided through by a 
normalized Euclidean model (i.e., [N mag' -Meg-^/lO - 6 '™- 12 - 5 )). Solid lines show an error-weighted least-squares fit to the counts. The 
calculated slopes are listed in the plots. 



normalized Euclidean model. For this global average of 
galaxy counts in J, H , and K s we find a = 0.59, 0.59, and 
0.63 (~ ±0.01 in all three bands), respectively. These 
values are all 0.01-0.02 above the Euclidean prediction, 
representing a possible galaxy space density increase to- 
ward fainter magnitudes of 7-15% over the magnitude 
range sampled. 

Next, we extracted 2MASS counts for areas of ~ 25 
deg 2 centered on each of our fields. In Figure [5] we show 
our raw counts (red diamonds) for JHK S < 17.5 plot- 
ted with 2MASS counts (black asterisks) from each of 
these 25 deg 2 subfields. The dashed lines show an error- 
weighted least-squares fit to the data. In the CLANS 
and CLASXS fields we include our counts extracted from 
the 2MASS-6x survey (green triangles), because these 
two fields lie within the 2MASS-6x Lockman Hole. The 
counts have been divided through by a normalized Eu- 
clidean model (a = 0.6). The calculated slopes are tabu- 
lated in the plots. We note that in many cases fitting to 
only the data points from this study would yield a much 
steeper slope than that obtained with the inclusion of 
2MASS. In this measurement of the local slope at our 
field positions, we find a mix of values ranging from sub- 
Euclidean to super-Euclidean. We note that the most 
sub-Euclidean slopes are found in our fields near the su- 
pergalactic equator (CLANS and CLASXS). We return 
to this point in Section |6~2"1 



In Figure [9l we show the error- weighted average of our 
counts (red squares) alongside the error- weighted average 
of the 2MASS counts (black diamonds) extracted from 
the 25 deg 2 subfields centered on our fields (for a total 
of 125 deg 2 ), as well as the counts from the 2MASS- 
6x Lockman Hole (green triangles). Again, we fit the 
slope with an error-weighted least-squares method and 
the slopes are listed in the plots. We find that, on av- 
erage, when combined with 2MASS, our counts yield 
a Euclidean slope within our error bars of ±0.01. We 
note that taken on their own, our counts in both the 
H and K s — bands would appear super-Euclidean. This 
is likely due simply to a lack of bright objects in the 
fields (HK S ~ 15), leading to large random errors in the 
brightest bin. 

6.2. Galaxy Counts and the Supergalactic Plane 

We now consider the implications of our bright galaxy 
counts in terms of local large scale structure. The 
most striking feature in the local large scale structure 
is the supergalactic Plane discovered bv lde Vaucouleursl 
[1951 The plane is roughly perpendicular to our Galac- 
tic plane. In Figure [SJ we find that in the CLANS 
and CLASXS fields, the slope appears distinctly sub- 
Euclidean, whereas in the other three fields the slopes 
are closer to the Euclidean prediction. This is interest- 
ing because the CLANS and CLASXS fields are of nearly 
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Fig. 8. — (a-e) J— band counts for 2MASS (black asterisks) vs. counts from this study in each field (red diamonds). Green triangles 
represent the counts from the 2MASS-6x survey, which overlaps with the CLANS and CLASXS fields, (f-j) J?— band counts in the same 
format, (k-o) K s — band counts in the same format. The blue dashed lines represent an error- weighted least squares fit to the data, and 
the calculated slope is listed in each plot. The counts have been divided through by a normalized Euclidean model to expand the ordinate. 



equatorial supergalactic latitude (—1.8 and —3.9 deg, re- 
spectively; see Table [T]). In fact, all of our fields are 
at relatively low absolute supergalactic latitude (average 
|SGB| ~ 10 deg), with A370 being the furthest off the 

plane at —25. 7 deg. 

The fields of lHuang et al.l (|1997f) . where a strong super- 
Euclidean slope was found in the K s — band, span abso- 
lute supergalactic latitudes from 10 < SGB < 54 with 
an average o f SGB ~ 30 deg. As a verification of the 
iHuang et all (|1997t ) result, we did another comparison 
in which we extracted 2MASS galax i es fro m patches of 
25 deg 2 centered on the IHuang et alj (|1997l ) fields. First, 
however, we extracted 2MASS counts from the actual 
areas of the IHuang et al.l (|1997| ) fields and found a ze- 
ropoint offset of —0.3 to their counts. We show this in 
Figure ITOl We applied this offset to their counts before 
proceeding with the comparison to the larger (25 deg 2 ) 
fields. 



We plot the zeropoint-adjusted counts from 
IHuang et all (|1997ft with the average of 2MASS 
counts over the larger surrounding areas in Figure 1111 
We calculate a slope of a = 0.67 for the combined 
dataset. Thi s value is slightly s maller than that found 
using just the lHuang et alj (|1997l ) datapoints (a = 0.69), 
but still distinctly super-Euclidean. Over the magnitude 
range sampled, this would still represent roughly a 
factor of two increase i n the space density of galaxies as 
f ound bv IHuang et al.l (|1997l ). 

iGardner et al.l (|1993f l found a super-Euclidean bright 
counts slope of a — 0.67 in their survey of a collection 
of fields t otalling ~ 10 de g 2 at a n average absolute SGB 
of ~ 25. Vaisanen et al. (2000) found a sub-Euclidean 
slope at bright magnitudes in the K s — band at high su- 
pergalactic latitude (SGB = 40), but this can likely be 
explained by the fact that they had three Abell clusters 
in their fields (A2168 at z = 0.06; A2211 at z = 0.15; 
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Fig. 9. — Average J, H, and if s -band counts for 2MASS sub- 
fields of 25 deg 2 each centered on our fields (black diamonds), 
2MASS-6x Lockman Hole (green triangles), and average counts for 
this study (red squares). The counts have been divided through 
by a normalized Euclidean model to expand the ordinate. 
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and A2213 at z = 0.15), which likely dominate their 
galaxy counts at bright magnitudes. The effect of just 
one galaxy cluster in a similar sized field can be observed 
in our data in Figure[TJ where the A370 (z = 0.37) counts 
show a large excess in the range 17 < JHK S < 19. The 
three lower redshift clusters in the lVaisanen et aTl (2000) 
fields would produce large excesses at brighter magni- 
tudes, and indeed they find their counts to be in excess 
over all other published studies at bright magnitudes. 

Thus, from the evidence above, it seems that a large 
fraction of the variability in the slope and absolute num- 
bers between different bright counts studies can be at- 
tributed to nearby large scale structure, be it individual 
galaxy clusters or local superstructure. When individ- 
ual nearby galaxy clusters are avoided, there appears to 
be evidence for the slope of bright galaxy counts to be 
related to supergalactic latitude. 

To investigate this further, we calculated the bright 
counts slope in 2MASS as a function of position on the 
sky by dividing the |b| > 30 catalog into slices in super- 
galactic latitude ranging in area from ~ 500 — 2500 deg 2 
each. As a result, we obtained slope measurements in 
J, H, and K s for -60 < SGB < 60. We plot the aver- 
age 2MASS slope as a function of SGB in Figure fT2l In 
this figure the Euclidean slope prediction has been sub- 
tracted from the measured values such that the dashed 
horizontal line represents the expectation for a Euclidean 
slope in all three bands. J-band is represented in blue, 
H— band in black, and K s — band in red. The error in 
the slope measurements is ~ 0.01 in all cases. We note 
that the similarity in slope measurements between bands 
suggests that our slope determinations are robust. 

In Figure [T^] the only distinctly sub-Euclidean slopes 
are found near the supergalactic equator, consistent with 
results from our CLANS and CLASXS fields, with all 
other latitudes appearing to be super-Euclidean. At high 



Fig. 10. — (a) The Hua ng et al.l U997I) K s counts are shown 
as blue diamonds. The 2MASS counts were constructed for the 
IHuang et all l(l997h fields for comparison and are shown as red 
asterisks. After applying azeropoint correction (K s = K B — 0.3; see 
text) the Huang et al. counts (black diamonds) align well with the 
2MASS K s counts from the same area. We adopt this adjustment 
as a correction to the Huang et al. counts, bringing them into 
good agreement with other studies over the same magnitude range. 
Error bars show 1 <r Poisson fluctuations, (b) The same data as 
in (a) but divided through by a normalized Euclidean model (i.e. 
[Nmag-Meg-^/lO - 6 ' 14 "- 12 - 5 )). 

northern supergalactic latitudes we find that slopes ex- 
ceeding t he Euclidean predict ion b y > 0.05, such as thos e 
found bv IHuang et al.l (|1997l ) and lGardner et all (|lQ93f ) . 
are typical. One interpretation of these results would 
be that near the supergalactic plane, the very brightest 
counts are populated with an excess of sources causing 
the slope to drop, and above the plane a relative paucity 
of nearby sources above the plane causing the slope to 
steepen. The southern supergalactic latitudes do not dis- 
play the same sharp rise to super-Euclidean as in the 
north. The magnitude range sampled in these counts is 
10.5 < JHK S < 15.5, which means that the brightest 
magnitude bin is centered on typical galaxies at a dis- 
tance of ~ 50 Mpc. This means that toward high north- 
ern SGB, where slopes are ~ 0.1 above the Euclidean 
prediction, the space density of galaxies at a distance of 
50 Mpc could be low by roughly a factor of two com- 
pared to the space density 250 — 350 Mpc distant. In the 
southern supergalactic cap the slopes are not as steep 
(~ 0.03 above Euclidean), but still represent a possible 
underdensity of ~ 25% at 50 Mpc relative to 250 — 350 
Mpc. 

Next we took a broad average of 2MASS counts in 
and out of the supergalactic plane. In Table [7] we give 
the measured average slope for galaxy counts in the su- 
pergalactic plane (|SGB| < 10) and out of the plane. 
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Fig. 11.— (a) The Huang ct al. (1997) K s counts are shown as 
black diamonds after a zeropoint correction (K s — 0.3; see text). 
The 2MASS counts were constructed from 25 deg 2 subfields (375 
deg 2 total) centered on the Huang et al. (1997) fields for compar- 
ison and are shown as red asterisks. The line shows an error- 
weighted least-squares fit to the data. We find a = 0.67 in this fit , 
a slightly lower number than was found for the [Huang ct al] (119971 ) 
data alone (a = 0.69), but still distinctly super-Euclidean. 



Fig. 13.— The 2MASS |b| > 30 sky was divided up into 
increasingly larger patches of area on the sky ranging from ~ 
25 — 10000 deg 2 to test for when the rms variability in slope be- 
tween sub-areas was reduced to the order of the variability induced 
by the supergalactic plane alone. The rms variability in measured 
slope is shown as a function of the area of the patches of sky. The 
result is that for patches of sky greater than ~ 1000 deg 2 , the 
variability in the slope can be accounted for by the supergalactic 
plane alone, such that other local large scale structure has been 
sufficiently averaged out. 



0.20 




-o.o5 r ... i ... i ... i ... i ... i ... i 

-60 -40 -20 20 40 60 

Supergalactic Latitude 



Fig. 12. — Average slope as a function of supergalactic latitude 
(SGB) for the 2MASS counts (|b| > 30; 10.5 < JHK S < 15.5). 
The expected Euclidean slope has been subtracted from the mea- 
sured values such that the dashed horizontal line at represents 
the Euclidean expectation in all bands. J— band appears in blue, 
H in black, and K s in red. Northern SGB slopes show a strong 
tren d toward super-Euc lidea n, generally consiste nt with results 
from lHuang et alj (IT997T) and lGardner et alj 1(1993?) . The only dis- 
tinctly sub-Euclidean slopes appear at SGB ~ 0, consistent with 
our findings in the CLANS and CLASXS fields which lie near the 
supergalactic equator. 

We find that on average the counts in the plane are 
sub-Euclidean, consistent with what we found in Fig- 
ure [TJ] The slope out of the plane is steeper by some 
Aa ~ 0.02 — 0.06 in both the northern and southern su- 
pergalactic caps, and super-Euclidean on average. Thus, 
the expected variability in the 2MASS bright counts 



TABLE 7 

2MASS Galaxy Counts Slope Over the Galactic 
Caps, the Supergalactic Equator and the 
Supergalactic Caps 



Region J Slope a H Slope a K s slope a 



Gal. |b| > 30 deg 0.59 0.59 0.63 

SGB| < 10 deg 0.55 0.57 0.59 

SGB< -10 deg 0.60 0.59 0.63 

SGB> 10 deg 0.61 0.61 0.65 

|SGB| > 10 deg 0.61 0.60 0.64 



a Errors in fitted slopes are ~ 0.01 in all cases 

slope is ~ ±0.04 due only to the very largest local struc- 
ture. 

We next divided the 2MASS |b| > 30 sky into 844 
patches of ~ 25 deg 2 each. We measured the slope of the 
galaxy counts from 11 < JHK S < 15 for each and found 
an rms variability in the slope of ~ 0.1. Inside distances 
of ~ 250 Mpc, 2MASS is 1 - 2 magnitudes deeper than 
M* and, as such, is measuring a fairly complete sample 
of the local luminosity function. This implies our 25 deg 2 
patches are sampling structure on size scales of ~ 20 Mpc 
and find rather strong variability. 

To test for the angular size over which all local large 
scale apart from the supergalactic plane could be suffi- 
ciently averaged out, we divided up the 2MASS |b| > 30 
sky into increasingly larger subfields and measured the 
rms variability in the measured slope across all subfields. 
Figure [13] shows the result of this exercise for subfields 
ranging from ~ 25 — 10, 000 deg 2 . We find that when av- 
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Fig. 14. — The slope of the 2MASS counts (black diamonds) in the range 11.5 < JHK a < 15.5 compared with the slope of the UKIDSS 
counts (red asterisks) in the range 13.5 < JHK S < 17.5. (a-c) show J— band counts, (d-f) show H— band counts, and (g-i) show K s — band. 
The top row of plots is for subfield 1, the middle for subfield 2, and the bottom for subfield 3 (see Table [8] for subfield details). The dashed 
and solid lines show an error-weighted least-squares fit to the 2MASS and UKIDSS counts, respectively. The slopes for each are given 
on the plots and fitted errors on the slopes are ~ 0.01 in all cases. The result is that for three large fields totalling ~ 135 deg 2 at high 
supergalactic latitude, no abrupt steepening of the counts slope is seen. 



eraging over subfields larger than ~ 1000 deg 2 , the rms 
variability between all subfields is of the same order of 
that induced by the supergalactic plane alone, and as 
such, all other local large scale structure has been suf- 
ficiently averaged out. In the 2MASS sample, this sug- 
gests that effects of local large scale structure have been 
sufficiently averaged out over scales of ~ 150 Mpc. 

Based on the analysis above, we conclude that the local 
supergalactic plane is overdense by 10 — 15% compared 
to regions at distances of ~ 250 — 350 Mpc, while the 
voids just above and below the plane are underdense by 
up to a factor of two relative to similarly distant regions. 

Co ld Dark Matter sim ulations such as the Millenium 
Run (Springcl et al.ll200"5l ) predict that the largest struc- 
tures in the universe should be < 100 Mpc in extent. 
Observed large scale structures, most notably the Sloan 



Great Wall (jGott et al.ll2005l ). demonstrate the existence 
of structure on scales much larger than simulations pre- 
dict. Consensus has not yet been reached on whether 
a structure like the Sloan Great Wall represents an ex- 
treme non-linearity in the matter distribution of the 
universe, or if inhomogeneities on several hundred Mpc 
scales are typical. Two recent analyses of the distribu- 
tion of matter in the S PSS DR6 dataset (jSarkar et al.l 
l2009t iKitaura et al.ll2009l ) demonstrate the current range 
of inte rpretation of observations, in which Sarka r~et al.l 
(120091) arrive at sca le of homogeneity of 70 Mpc and 
IKitaura et alJ <|2009T l claim to have discovered a new large 
void 150 Mpc in diameter. Clearly the typical size of 
large scale structure still depends on how one interprets 
the data. In our analysis, we find that averaging galaxy 
counts over areas on the sky spanning ~ 150 Mpc suf- 
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TABLE 8 

Three Subfields Selected from the UKIDSS DR3 LAS 



Subfield 


ARA (deg) 


ADEC (deg) 


ASGB (deg) 


1 


201.5 - 222.1 


9.0 - 11.4 


10 - 30 


2 


337.5 - 360.0 


-1.0 - 1.0 


13 - 34 


3 


135.0 - 150.0 


7.0 - 11.4 


-(37 - 52) 



ficiently smooths out local structure. However, in the 
same data, the super-Euclidean slopes above and below 
the supergalactic plane suggest that the local space den- 
sity of galaxies may be low by 25 — 100% compared to 
regions a few hundred Mpc away, which would suggest 
that local structure exists on much larger scales than 150 
Mpc. 

6.3. Comparison with the UKIDSS Large Area Survey 

We employ the UKIDSS survey to investigate the pos- 
sibility that the entire 2MASS volume resides in a rel- 
atively underde nse region th a t is > 250 Mpc in radius, 
as suggested bv IHuang et al.l (|1997l ). The expectation if 
this were the case would be a steepening of the galaxy 
counts slop e in the range 15 < J HK S < 17, as, in fact, is 
seen in the IHuang et alj (|1997l ) sample when compared 
with 2MASS. 

As a test for a lar ge void, we employ the UKIDSS 
(jLawrence et al.l [20071 ) third data release (DR3; Warren 
et al., in preparation) Large Area Survey (LAS), which 
has recently become public. The DR3 LAS covers sev- 
eral hundred degrees of sky within the SDSS footprint to 
depths of JHK S ~ 19. We selected three large subfields 
totalling - 135 deg 2 in the UKIDSS LAS. We down- 
loaded catalogs from the WFCam Science Archive, and 
the range in celestial coordinates and SGB for these sub- 
fields are given in Table [8] We selected the fields to span 
a wide range of SGB above and below the supergalactic 
plane. 

We cross-correlated the UKIDSS catalogs with their 
SDSS counterparts and removed stars from the catalogs 
via the SDSS classifier and g'JK color-color diagrams in 
a similar manner to what is described in Section 15.11 A 
comparison between the slope of the 2MASS counts and 
the UKIDSS counts for the three subfields is shown in 
Figure [TU 2MASS counts are shown in black diamonds 
and UKIDSS counts in red asterisks. The dashed and 
solid lines show error-weighted least-squares fits to the 
2MASS and UKIDSS counts, respectively. The slopes 
for each fit are listed in the plots, and the associated 
errors in the fitted slopes are ~ ±0.01 in all cases. 

The result is that for these three subfields, we do not 
observe a steepening of the slope at magnitudes beyond 
the 2MASS limits. In fact, in the majority of cases in 
Figure [T31 the slope appears to be decreasing into the 
fainter magnitudes of the UKIDSS survey. 

7. SUMMARY 

We have presented a deep, wide-field NIR survey over 
five widely separated fields at high galactic latitude cov- 
ering a total of ~ 3 deg 2 in J, H, and K s . The deepest 
parts of the survey (~ 0.25 square degrees) reach 5 a 
limits JHK S > 24. As such, this is one of the deepest 
wide-field NIR imaging surveys to date. In this paper we 
focus on the bright galaxy counts slope and the implica- 



tions for local large scale structure. We leave analysis of 
the faint galaxy counts and other applications of these 
data for future papers. 

We measure the slope of the bright galaxy counts in 
our data combined with the larger 2MASS fields at our 
positions on the sky, and we find it consistent with the 
Euclidean prediction, on average, over our five fields. We 
note that our fields are at a low average supergalactic 
latitude (SGB ~ 10) and our fields near the supergalac- 
tic equator show a sub-Euclidean slope. In the 2MASS 
|b| > 30 sky, we find the slope of the counts to be 
sub-Euclidean near the supergalactic equator and super- 
Euclidean at all other SGB. This is consistent both with 
our measured slope near the supergalactic equator and 
with other studies that have found steeper slopes above 
the plane, except when nearby galaxy clusters are in- 
cluded in the counts. This can be understood, at least in 
part, in terms of the fact that our galaxy counts begin- 
ning at 11 th magnitude are sampling the structure of the 
plane itself in sightlincs along the plane, and sampling 
the voids above and below the plane in sightlincs away 
from the plane. 

We further explore local large scale structure in the 
2MASS sample by varying the area over which the counts 
are averaged to find the angular scale on which the vari- 
ability between one area and another is of the same order 
as that expected to be induced by the supergalactic plane 
alone. We find that fields of > 1000 deg 2 achieve this re- 
sult, corresponding to averaging over scales of ~ 150 Mpc 
in the 2MASS volume. This is to say, we find that on 
scales of ~ 150 Mpc, local large scale structure is suffi- 
ciently averaged out of the counts. However, our result 
that the galaxy counts slope is super-Euclidean above 
and below the supergalactic plane implies that the lo- 
cal space density of galaxies (away from the plane itself) 
could be low by 25 — 100% relative to regions a few hun- 
dred Mpc distant. This suggests that local structure may 
exist on scales much larger than 150 Mpc. 

Finally, we explore the possibility of whether the entire 
2MASS volume exists in a region of relative underdensity 
through a comparison with the UKIDSS DR3 LAS. We 
find that the slope of the UKIDSS counts is either con- 
sistent with 2MASS or takes a lower value over the three 
large subfields we selected. This result is inconsistent 
with the expectation of a steepening of the counts slope 
in the range 15 < ,JHK S < 17 that would be expected if 
we lived inside a large void of radius ~ 300 Mpc. 

Above all else, the comparison and analysis of several 
large NIR surveys presented here has served to highlight 
the complexity of deciphering local large scale structure 
from galaxy counts alone. The completion of surveys 
like UKIDSS will expand our view of the local Universe, 
but as demonstrated here, even surveys of large areas 
on the sky may suffer biases due to local structure. As 
such, mapping out local large scale structure remains a 
complex problem that will likely remain ambiguous until 
extensive spectroscopy in combination with large scale 
photometric surveys can generate a three dimensional 
picture of the local Universe. 
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